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A VARIATIONAL PRINCIPLE FOR MAGNETOHYDRODYNAMIC CHANNEL FLOW 

by Norman C. Wenger 
Lewis Research Center 

SUMMARY 

A variational formulation is presented for a class of magnetohydrodynamic (MHD) 
channel flow problems. This formulation yields solutions for the fluid velocity and the 
induced electric potential in a channel with a uniform transverse static magnetic field. 
The channel cross section is constant but arbitrary, and the channel walls can be either 
insulators or conductors with finite electrical conductivity. Electric currents are per- 
mitted to enter and leave the channel walls so that the solutions are suitable for MHD 
generator and pump applications. An example of a square channel with conducting walls 
is solved as an illustration. 


INTRODUCTION 

The study of magnetohydrodynamic (MHD) channel flow has received considerable 
attention in the past decade. This interest. has been motivated by three principle applica- 
tions; the MHD generator, the MHD pump, and the electromagnetic flowmeter. 

The general model that is normally considered in these studies consists of an infi- 
nitely long channel of constant cross section with a uniform static magnetic field applied 
transverse to the axis of the channel. The walls of the channel are either insulators, 
conductors, or a combination of insulators and conductors depending on the intended ap- 
plication. 

For example, in the MHD generator and pump cases, the channel cross section is 
normally rectangular with insulated walls perpendicular to the magnetic field and con- 
ducting walls parallel to the magnetic field. For the electromagnetic flowmeter case, 
the channel cross section is normally circular with conducting walls. 

In order to carry out an analytical solution for MHD channel flow, it is generally 
necessary to make simplifying assumptions such as requiring the channel walls to be 
either perfect conductors or perfect insulators or requiring the channel walls to be very 
thin. These and other simplifications often greatly limit the usefulness of the results. 



This is particularly true for the electromagnetic flowmeter case since the thin wall ap- 
proximation is often not valid for liquid metal applications and the wall conductivity is 
neither zero nor infinite. In addition, many analytical solutions give results in the form 
of infinite series which converge poorly for the large values of the static magnetic field 
that are encountered in practice. 

To alleviate some of these difficulties, Tani (ref. 1) developed a variational formu- 
lation for the solution of MHD channel flow problems. His formulation gives solutions 
for the velocity profile and the induced magnetic field distribution in the channel for an 
arbitrary channel cross section. It requires, however, that the channel walls be either 
perfect conductors or insulators and that the admissible functions for the velocity and 
induced magnetic field satisfy appropriate boundary conditions. 

In this report, a variational formulation is presented that gives solutions for the 
velocity profile and the electric potential distribution in a chanhel of arbitrary cross sec- 
tion. It also gives solutions for the electric potential distribution in the channel walls. 
The walls of the channel can be a combination of insulators and conductors but the con- 
ductors may have a finite conductivity. In addition, the admissible functions for the 
velocity and potential need not satisfy any prescribed boundary conditions. Moreover, 
the formulation is sufficiently general to allow electric currents to enter and leave the 
channel walls so that the solutions obtained are suitable for the MHD generator and pump 
applications. 

The report concludes with an example that consists of a square channel with conduct- 
ing walls of finite conductivity. 


THE MODEL 

A cross section of a generalized channel is shown in figure 1. It consists of the 
fluid duct bounded by the conducting walls S c and the insulated walls The con- 
tours C^ c and denote the fluid- conducting wall interface and the fluid- insulated wall 
interface, respectively. The contour C CQ denotes the outer edge of the conducting wall. 
The vector n is the unit normal to the contours with the positive direction as shown. 

The applied static magnetic field B Q is parallel to the x-axis and is uniform with 
respect to y and z. The applied or generated current density at the outer edge of the 
conducting wall J is considered positive when directed outward. It is assumed that the 
net current entering the channel cross section due to J & is zero, so that the two- 
dimensional features of the model are retained. 
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BASIC EQUATIONS 

The basic equations to be used are the standard MHD equations for steady- state 
conditions which consist of Maxwell's equations, the momentum transport equation, and 
the generalized Ohm’s law. These are 


V x E = 0 

(la) 

V x 5 = p Q J 

(lb) 

v . 5 = o 

(1c) 

p (v . v)v = -v p + j x i + 77 V 2 v 

(Id) 

J = a f (E + V X 5) 

(le) 


where E, B, J, and p Q are the electric field intensity, magnetic flux density, electric 
current density, and magnetic permeability of free space, respectively; and V, p, rj, 
cr f , and p are the fluid velocity, density, viscosity, electrical conductivity, and pres- 
sure, respectively. 

Equations (la) to (le) are based on the assumptions that the fluid is homogeneous and 
incompressible, the magnetic permeability of the fluid is the same as that of free space, 
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the electric charge density and Hall current are negligible, and the fluid flow is laminar. 
It may be argued that the last assumption of nonturbulent flow greatly limits the useful- 
ness of the results since many flows in practice are turbulent. However, it has been 
found experimentally (ref. 2 ) that the onset of turbulence occurs at a much higher 
Reynolds number in MHD flow than for ordinary flow due to the suppression of turbulence 
by the static magnetic field. 

The five basic equations (la) to (le) can be combined to give two coupled equations of 
second order. First, an electric potential U which satisfies equation (la) identically 
can be defined as 


E = -VU (2a) 

Next, substituting J from equation (le) into (lb) and then taking the divergence of equa- 
tion (lb) give 


V 2 U = V • (V x 5) (2b) 

where equation (2a) has been used to eliminate E. The second equation can be obtained 
by substituting J from equation (le) into (Id) and then eliminating E using equation (2a) 
giving 


p(V . V)V = -Vp - cr f [vu X B - (V X 5) X 5] 


? 7^ 2 V 


( 2 c) 


Due to the uniformity of the channel cross section and the applied magnetic field with 
respect to the z-axis, all quantities in the basic equations are independent of z except for 
the pressure which is linear in z (ref. 3). In addition, it can be shown that the fluid ve- 
locity V has only a z- component V and that E and J have only x- and y- components 
(ref. 4). Furthermore, the total magnetic field B consists of the applied field B Q in the 
x-direction and an induced field B^ in the z-direction. 

Since the velocity does not vary with z and has only a z- component, the first term in 
equation (2c) vanishes. Expanding the vector cross products and taking the z- component 
of equation ( 2 c) along with equation ( 2 b) give the following governing equations: 


S 2 U t 9 2 U b 

ax 2 ay 2 " • ° 


(3a) 


a 2 v„ 


ax' 


a 2 v 


sy 


i 

r 7 3z 


3p , g f B 0 3U 


v 3y 


V 


v z = 0 


(3b) 
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These equations apply, of course, only in thb fluid duct region Sj. In the conducting wall 
region S c , equation (3a) applies with V z = 0. In the insulated wall region S^, equa- 
tion (3a) also applies with V = 0 but it need not be solved. 

Equations (3a) and (3b) are not unique governing equations in the sense that variables 
other than the velocity and electric potential can be selected for retention. Tani (ref. 1), 
for example, eliminated the electric potential but retained the velocity and the induced 
magnetic field yielding a different but equivalent set of governing equations. 

In addition to the basic equations, appropriate boundary conditions must be specified 


to determine the solution uniquely. These conditions are 

V z = 0 on C fc and C f . (4a) 

f 

U - U = 0 on C fc (4b) 

f w 

OjVU • n - a w VU • n | = 0 on C^ c (4c) 

f w 

a f vu • n = 0 on C fi (4d) 

f 

CT W VU ‘ » + J a = 0 on C co (4e) 

w 


where | and | refer to evaluating the quantity on the fluid or wall side of the contour, 

f w 

respectively. 

The boundary condition equations (4a) to (4e) require the following: 

(1) The fluid velocity must vanish on the fluid-wall interfaces C^ c and C^. 

(2) The electric potential must be continuous across the fluid- conducting wall inter- 
face . 

(3) The component of the electric current normal to the fluid- conducting wall inter- 
face Cj c must be continuous. 

(4) The component of the electric current normal to the fluid- insulated wall inter- 
face must vanish. 

(5) The component of the electric current normal to the outer edge of the conducting 

wall C! must equal the applied or generated current J 
co a 

In solving the equations, it is convenient to work with dimensionless quantities. This 
can easily be accomplished by defining L and V Q to be a characteristic length and 
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characteristic velocity of the channel. Let 


X = — , Y = Z =2. 


dimensionless coordinates 


B LV 
o o 


dimensionless potential 


V = — — 


dimensionless velocity 


M = B q L 


Hartmann number 


P - ~L 9p 

° W Q 3z 


dimensionless pressure gradient 


B V a 
o o w 


dimensionless applied or generated current 


ratio of wall-to-fluid conductivity 


Combining equations (5a) to (5g) with equations (3a) and (3b) yields the following set of 
equations in dimensionless form: 


9 2 W f 9 2 W 9V _ 0 
ax 2 9Y 2 9Y 


on S r 


i!v + 3?V +p 2 3 W_ m 2 v = 0 

9X 2 9Y 2 0Y 


on 


9 2 W + 9^W = q 
9X 2 9Y 2 


on S, 


I II II ■ II II I I 


Likewise, combining equations (5a) to (5g) with equations (4a) to (4e) gives the following 
set of dimensionless boundary condition equations: 


V | 

f 


= 0 


on C fc and C fi 


(6d) 


W 



= 0 


on C fc 


(6e) 


i /dY aw aw\ 

y 

/ . \dx ax ay/ 

V-(gJ 

,1 /dy\^ 

f i/ 1 +(— ] 

T \dX/ 


dY aw 
.dx ax 



on C fc (6f) 




dY aw 
dx ax 



- J o = ° 


w 


on 


co 


(6g) 


(6h) 


The unit normal vector 


n has been replaced by 


n = 


dY 

dX 


a x + a y 



where a x and Sy are the unit vectors in the X- and Y-directions, respectively. The 
sign of the square root must be selected so that the positive direction for n is as shown 
in figure 1. 


VARIATIONAL EXPRESSION 

The goal of this section is to construct a functional of the dependent variables V 
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and W so that the associated Euler-Lagrange equations are the basic governing equations 
(6a) to (6c) and where the corresponding natural boundary conditions are the prescribed 
boundary condition equations (6d) to (6h). This construction is performed by summing 
terms that are obtained by multiplying each governing equation and boundary condition 
equation by a suitable function and then integrating over the corresponding area or contour 
where the equation is valid. 

Let 6V and 6W be arbitrary functions of X and Y that are continuous with piece- 
wise continuous first derivatives. The integrals 



are identically zero for any 6V and 6W since the quantities in brackets [ ] are identi- 
cally zero by virtue of equations (6a) to (6c). 

Four additional integral expressions that are identically zero can be obtained from 
the boundary condition equations (6d) to (6h) by integrating along appropriate contours. 

Recalling that the differential length along a contour is given by y 1 + (dY /dX )^ dX, these 
integral expressions can be defined as 
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These integrals vanish since the quantities in brackets are zero by virtue of the boundary 
condition equations (6d) to (6h). 

The symbol 6 can be defined as the variation operator so that the functions 6V 
and SW can be considered as the variation of V and W, respectively. In addition, the 
6 operator commutes with 0/9X and 9/9Y since X and Y are independent variables. 
Using these properties along with Green's lemma, the integrals 1^ to 1^ can be inte- 
grated by parts and combined to give 


7 


£ 


6F 


where 


(8a) 



Since each I (n=l, . . 7) is identically zero, the variation of the functional F is 

zero. Thus, F is stationary; that is, first order changes in V and W about their true 
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values produce only second order changes in F. 

Another way of stating this result is that of all functions that are continuous with 
piecewise continuous first derivatives, the particular pair of functions for V and W 
that make F stationary, satisfy both the basic equations (6a) to (6c) and the boundary 
condition equations (6d) to (6h) and, hence, are the desired solutions. 

Even though F was shown to be stationary, it does not necessarily mean that F 
has a maximum or minimum at the true solution for V and W. For example, as V 
and W are varied from their true values, F may always increase, always decrease, or 
either increase or decrease depending on how V and W are varied. To determine 
which case corresponds to the F under consideration the quantity F(V + 5V, W + 6W) - 
F(V, W) is computed giving 


F(V+ 6V, W + 5W) - F(V, W) = 2 



- m 2 v + m 2 £5^ 

(Lax 2 ay 2 

3Y 


6V 
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A careful examination of equation (9) reveals that the first six integrals vanish be- 
cause the quantities in brackets [ ] are identically zero. The remaining integrals are of 
second order in 6V and 6W. This result is not surprising since F was constructed so 
that all first order terms in SV and SW vanished. Of the three remaining integrals in 
equation (9), two are negative definite and the last can be of either sign. Thus, F has 
neither a minimum nor maximum at the true values for V and W. However, if the class 
of admissible functions for V and W is restricted so that the last integral must vanish, 
then F corresponds to a maximum at the true values for V and W since then 
F(V + 6V, W + 6W) - F(V, W) < 0. 

A study of the last integral in equation (9) reveals that the proper restriction to im- 
pose is that V must vanish on the contours and C^. An alternate choice which 

also makes the last integral vanish is to specify 3V/Sn on C^ c and C^. This choice is 
useless, however, since it would require solving the problem another way first to deter- 
mine the correct value for 3V/9n. 

Requiring V to vanish on C£ c and may provide a great simplification in many 
problems in obtaining approximate values for V and W since finding a maximum for F 
is often much easier than finding a stationary point. Moreover, requiring V to vanish on 
Cj c and completely eliminates one integral in the expression for F given by equa- 
tion (8b). 

Before a solution for V and W can be determined, values for P Q and J Q must be 
specified. The dimensionless pressure gradient P Q must be a constant as previously 
noted. The dimensionless current density J Q , however, can be specified as a function of 
the coordinates along the contour C cq . Since the basic equations and boundary condition 
equations are linear in V, W, P Q , and J , solutions for V and W can be obtained by 
superimposing solutions for P Q ^ 0 and J Q = 0 with those for P Q = 0 and J Q 4 0. 

To complete the study of the variational expression, it is desirable to determine its 
physical significance. 


PHYSICAL SIGNIFICANCE OF VARIATIONAL EXPRESSION 


Consider the power or energy balance that exists in MHD channel flow. The power 
per unit length that is supplied to the fluid by the pressure gradient P^p is given by 


P 




dx dy = 7 )V° 



P Q V dX dY 


(10a) 


If this quantity is negative, it simply means that the channel is acting as a pump. 
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The dissipative terms consist of the viscous losses in the fluid P^, the ohmic losses 

in the fluid P , and the ohmic losses in the conducting walls P . Expressing each of 
°f °w 

these in terms of the power dissipated per unit length along the channel gives 


/ 


V # 1 


ra v_ a v, 

— + - 


L0X 3y 




21 V„ dx dy = ,v 2 / |(!i) + (^ 

3Y 


dX dY 


r?V 


/ 

^Cfc^fi 


V( IX dY - 9V dx 
3X 3Y 


/ 


V'*/ 


3x A 3y 


dx dy 


(10b) 


’I [« 


= , v > 2 / \m\\(™-v) 2 


3X / \3Y 


dX dY 


(10c) 


P = CT 
a w w 


/ 


?£S -™a 

ax x ay y 


dx dy 


2 2 

= i7V‘y»T 


/ 


aw\ 2 /aw x2 


ax / \3Y 


dx dY 


(lOd) 


Equation (10b) reveals that the viscous losses in the fluid can be split into two parts; 
the volume losses P and the surface losses P„ where 


V 


Vc 
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The surface losses are zero when the boundary condition is imposed that requires V to 
vanish on Cf and C^. 

The remaining term to be considered is the loss due to the current J . Since J is 
positive by definition when it is directed outward, the power that is supplied to an ex- 
ternal load per unit length of the channel Pj is given by 



If P T is negative, it simply indicates that power is being supplied by an external source. 
J a 

Since power is conserved, the power balance for the channel can be expressed as 


P = P 
Ap *V 


+ P - + P + P 


+ P, 


(ID 


w 


Comparing the expression for the functional F given by equation (8b) with the various 
power dissipation terms given by equations (10a) to (lOg) reveals that F can be ex- 
pressed as 


t?V F = 2P A - P„ - P - P„ 

0 A P a f CT w 


- 2P „. - 2 p j, 


( 12 ) 


A word of caution is in order at this point. The expression for F given by equa- 
tion (8b) is defined and valid for an arbitrary choice for V and W. Likewise, the power 

dissipation terms P , P , etc. , given by equations (10a) to (lOg) are valid for arbi- 
V v a f 

trary values for V and W. Thus, equation (12) is valid, in general. However, equa- 
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tion (11) which is the power balance for the channel is only valid for the correct solutions 
for V and W. 

As shown in the previous section, the stationary point for F corresponds to the true 
solutions for V and W. For these values only, equations (11) and (12) can be combined 
to give 


F 


st 



(13) 


where it has been recognized that P vanishes for the time V. Thus, the stationary 

”s 

value of F is proportional to the difference between the power supplied to the fluid by 
the pressure gradient and the electrical power delivered to an external load. 

An important special case occurs for J = 0 which yields 

a 


st 


-I 


P V dx dY 
o 


(14) 


where P^ has been replaced using equation (10a). Since the dimensionless pressure 
gradient P Q is a constant, the stationary value for F is proportional to the average 
fluid velocity in the channel. This is a very important result since the average velocity, 
which is often the main quantity of interest, is proportional to a stationary quantity which 
can be computed to good accuracy. 

If the boundary condition V = 0 on Cj c and is satisfied by all admissible 
functions, F has a maximum at its stationary point as noted in the previous section. The 
maximum for F using a subset of the class of admissible functions for V and W will 
be less than the maximum for F using the entire class of admissible functions. Thus, a 
lower limit for the average velocity can be easily found by using any admissible function. 


EXAMPLE: SQUARE CHANNEL WITH CONDUCTING WALLS 
Variational Solution 

A square channel is shown in figure 2 using the dimensionless coordinates. The 
characteristic length for the channel L has been chosen as one-half its height or width 
so that its inner corners are located at (1,1), (1,-1), (-1, 1), and (-1,-1). The normal- 
ized wall thickness t is the actual wall thickness divided by L. 
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1 + t 
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hi 


-l -t 


x 


Figure 2. - Cross section of square channel with conduct- 
ing walls. 


Approximate solutions for the velocity and electric potential will be obtained using 
the Ritz technique. In this technique, the velocity and potential are expressed in terms 
of known functions of X and Y that approximate the true solution but contain adjustable 
parameters Ap . . A . The approximate solutions for V and W are then substi- 
tuted into the expression for F given by equation (8b) and the indicated integrations with 
respect to X and Y are performed. This leaves F as a function of the parameters 
Ap . . . , A r and the characteristic parameters of the channel P , y, M, and J As- 
suming that V vanishes at X = ±1 and Y = ±1 for all values of Ap . . . , A n , the sta- 
tionary value of F can be found by maximizing F with respect to Ap . . ., A . The 
corresponding values for Ap . . . , A r when substituted into the approximate functions 
for V and W will yield the closest approximations to the velocity and potential that are 
possible for the class of functions used. 

In order to determine the accuracy of the solutions obtained using the Ritz technique, 
a sequence of approximations is normally used. In this method, a suitable, complete, in- 
finite set of functions is selected such as sine or cosine harmonics, Bessel functions, 
etc. , so that some linear combination of them is capable of representing the solution. 

For the MHD channel flow problem these functions must be capable of representing any 
continuous function with a piecewise continuous derivative over the channel cross section. 

The procedure consists of first obtaining an approximate solution using a linear com- 
bination of a finite number of these functions where the adjustable parameters Ap . . . , 
A n correspond to the coefficients of these functions. The problem is then repeatedly 
solved, each time increasing the number of functions used. By comparing the results 
from successive approximations, an estimate of the accuracy and convergence rate can be 
obtained. 
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In many problems, however, much is already known about the solution so that the 
trial functions can be tailored to more accurately approximate the solution thereby re- 
ducing the number of adjustable parameters required. This alternate procedure does not 
allow the error in the approximation to be easily estimated but it is much easier to use 
computationally since fewer adjustable parameters are involved. This alternate proce- 
dure will be followed in this example. 

The solution for the square channel will be determined for P = 1 and J = 0. Let 

oo 

the trial functions for V and W be given by 


V(X,Y) = Aj 




0 <X < 1 
0 < Y < 1 


W (X , Y ) 



+ c 2 Y 




0 < X < 1 + t 
0 £ Y < 1 + t 


(15a) 


(15b) 


where Ap Cp C 2 , Cg, ctp /3 2 , and are adjustable parameters. Because 

of the symmetry of the problem, it is only necessary to specify V and W in the first 
quadrant. For other quadrants, V and W can be found using the relations V(X,Y) = 
V(X, -Y) = V(-X,Y) and W(X,Y) = W(-X,Y) = -W(X,-Y). Since the admissible functions 
for V and W must be continuous with piecewise continuous derivatives, all exponents 
in equations (15a) and (15b) must be 1 or greater. 

Trial functions of the form given by equations (15a) and (15b) have proven to be very 
useful in solving MHD channel flow problems of this type since the velocity profile, for 
example, can go from a parabolic (o'j = o' 2 = 2)to nearly slug flow (aip a 2 large) by 
merely varying two parameters. 

Substituting equations (15a) and (15b) into the expression for F given by equa- 
tion (8b) yields after performing the integrations 

F=rf[A]W- [^] T [D] (16) 

where 



_A 1 

C 1 

[A] = 

a ll a 12 a 13 

a 12 a 22 a 23 

P>] = 

V 

0 


C 2 

L J 


_ a 13 a 23 a 33_ 


0 
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a ll= - 


4 


8 44 


(a 2 + 1)(2 Q' 1 - 1)(2q! 2 + 1) (a^ + l)(2a 2 - l)(2o' 1 + 1) 


16 M 2 a^Q !2 


(«! + l)(2a 1 + l)(a- 2 + 1)(2 af 2 + 1) 


a 12 = 


4 M 2 Q! 1 Q !2 

1 + 

C 3 

01 + Q -2 

+ 1 

(/ 3 g + l)(/ 3 g + Q!^ + 1 ) 


l 13 


4 M ^ 1 « 2 


H + a 2 


+ 1 (/ 3 g + + 1 ) 


a 22 * - 4M 


P 3 C 3 


4 


2 Co C, 
1 + d 


(2j3g - 1)(2^ +1) 2^-1 ^g + l 2^ + 1 


+ yt 


2 (/ 3 1 + | 3 3 )/ 3 |c| 2 Z 3 2 


( 2 / 3 3 - 1 )( 2 / 3 1 + 1 ) 2 / 3 j - 1 


; 0 i +- 


^ 1^3 + Vi + <V C 3 
0 3 + 1 2 / 3 3 + 1 


a 23 = “ 4M 


/ 3 2 C 2 

P 3^3 


0 n 0. 


+ --- 


1 F 2 


( 2 j 3 3 - lK/Sj + 0 2 + 1 ) fi 1 + / 3 2 - 1 


2 C q ci 
1 + — — + J 


/3 q + 1 2 / 3 ^ + 1 


+ yt 


(^1 + @2 + 2 ^ 3 ^ 3^3 


( 2 / 3 3 - 1 )(/ 3 1 H-fy+l) 


0 i 0 ^ 


1^2 


01 + 02 “ 1 


( 0 J + 0 £) + 


2(0j + 0 2 + 0 3 )C 3 (0 j + 0 2 + 2/3 3 )C 2 


0 3 + 1 


20 3 + 1 
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a 33 


= -4M 2 



0 2 C 2 

P 3 C 3 

- 1 )( 20 2 + 1 ) 




2C 3 
J3g + 1 



+ yt 


20 2 + /3 3 )/3|c| 
(2^ - l)(2j3 2 + 1) 




(2/3 2 + /3 3 )C 3 

V * 71 


(^2 + ^3)^ 

2 V^ r 



^ J 


d 


1 


8of l a 2 _ 

(«! + 1)(q , 2 + !) 


In obtaining equation (16), the "thin wall" approximation t « 1 was made so that 
the results are directly comparable with published values. This approximation did not 
have to be made in order to use this formulation but was done since only thin wall results 
are available for comparison. 

The maximum value for F and the corresponding values for the parameters were 
found using a computer. Since the normalized cross-sectional area of the channel is 4 
and P Q = 1, the proportionality constant between the average dimensionless fluid veloc- 
ity V and Fgj. is 1/4 (see eq. (14)). Values for V as a function of the Hartmann num- 
ber M are shown in figure 3 for various values of conduction parameter yt. 

The results show that the average fluid velocity decreases as M and yt increase. 
This decrease can be explained as follows: For M / 0, an electric current is induced in 
the fluid due to the fluid motion. If yt = 0, the wall is an insulator and, hence, any in- 
duced current must form a closed path entirely in the fluid. Thus, the total electromag- 
netic force on the fluid is zero. However, since the electromagnetic force is not iden- 
tically zero everywhere, the velocity profile of the fluid is distorted which, in turn, in- 
creases the viscous force and consequently reduces the average fluid velocity. For 
yt / 0, the induced current path is located partially in the wall so that there is a net elec- 
tromagnetic force in addition to the viscous force acting to reduce the fluid velocity. 


Comparison of Results 

The variational solution can be compared with other solutions for some limiting 
cases. For M - 0, the average dimensionless velocity from the variational solution is 
0. 1403, independent of yt, as compared with the exact value of 0.1406, which csm be 
computed using Fourier expansion techniques. 
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The exact solution for the average flow in a rectangular channel with insulated walls 
(yt = 0) and perfectly conducting walls (yt = °°) has been obtained by Shercliff (ref. 3) and 
Chang and Lundgren (ref. 5), respectively. Each of these solutions is in the form of a 
series which converges poorly for large M. Williams (ref. 6), however, transformed 
these solutions and obtained asymptotic forms for the average velocity for large M. 
These solutions, simplified for P Q = 1 and the square channel, are as follows: 



10' 1 10 ° 10 1 10 2 


Hartmann number, M 

Figure 4. - Comparison of variational solution with Williams (ref. 6) asymptotic solutions for average 
dimensionless velocity. 
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A comparison between the variational solutions and these asymptotic forms is shown 
in figure 4. The agreement between the two solutions for yt = 0 is excellent for M ^ 10. 
The variational solution for yt = 10 4 is always slightly less than the asymptotic form for 
yt = The difference, however, decreases to less than 0. 1 percent at M = 1000. 

A Fourier expansion type solution for the rectangular channel with thin walls of 
finite conductivity has recently been obtained by Chu (ref. 7). A comparison of his solu- 
tion with the variational solution is shown in figure 5. The agreement between these two 
solutions is also quite good. As shown, the variational solution for the average velocity 
is always slightly less than the series solution value. This is due to the fact that the 
computed maximum for F, and, hence, the average velocity, is always less than or equal 
to the true maximum for F since the trial functions used are a subset of the entire class 
of admissible functions. 



Figure 5. - Comparison of variational solution with Chu's 
(ref. 7) solution for average dimensionless velocity. 


21 



CONCLUDING REMARKS 


A variational formulation was presented for a class of MHD channel flow problems. 

A stationary expression was developed that yielded solutions for the fluid velocity and the 
induced electric potential in a generalized channel. An example of a square channel with 
conducting walls was solved as an illustration. Very good agreement was obtained be- 
tween the variational solution for the average velocity in the square channel and other 
published values. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 24, 1969, 

120-27. 
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APPENDIX - SYMBOLS 


[A] 

matrix (see eq. (16)) 

n 

unit normal vector to con- 

A 1 

adjustable parameter 


tours 

a ij 

A 

a , a 

ij L 1 element of matrix [A] 

G 

order of 

unit vectors in x- and 

p j a 

power delivered to ex- 

x’ y 

B 

y-directions 
magnetic flux density 

p 

A 

ternal load per unit 
length of channel 

dimensionless pressure 

B i 

induced magnetic flux 
density 

u 

P Ap 

gradient 

power supplied to fluid by 

B o 

uniform applied magnetic 
flux density 

pressure gradient per 
unit length of channel 

C fc’ C fi’ C co 

contours in generalized 
channel cross section 
(see fig. 1) 

P , 

power dissipated by vis- 
cous force per unit length 
of channel 

C 1’ C 2’ C 3 

adjustable parameters 

P u 

power dissipated by surface 

[D] 

matrix (see eq. (16)) 

'S 

viscous force per unit 

d l 

element of matrix [D] 


length of channel 

E 

electric field intensity 

V 

power dissipated by volume 
viscous force per unit 

F 

stationary functional 
(see eq. (8b)) 


length of channel 
power dissipated by ohmic 

F st 

stationary value of F 

a f 

loss in fluid per unit 

Ij, . . • > I rj 

integrals (see eqs. 


length of channel 


(7a) to (7g)) 

P n 

power dissipated by ohmic 

J 

electric current density 

w 

loss in walls per unit 

J a 

applied or generated elec- 


length of channel 


tric current density 

P 

pressure 

J o 

dimensionless applied or 
generated electric cur- 
rent density 

s f ,s , s. 

f’ c’ 1 

surfaces in generalized 
channel cross section 
(see fig. 1) 

L 

characteristic length of 

t 

dimensionless wall thick- 


channel 


ness of square channel 

M 

Hartmann number 


(see fig. 2) 
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u 

electric potential 

V 


fluid viscosity 

V 

dimensionless fluid 


' • X n 

generalized adjustable 

V 

velocity 

average dimensionless 

^o 

parameters 

magnetic permeability of 

V 

fluid velocity 
fluid velocity 


> 


free space 
fluid density 

V D 

characteristic velocity 

a f ,i 

T w 

fluid and conducting wall 

V z 

w 

z- component of V 
dimensionless potential 

y 

electrical conductivities 
a w/ a f 

X, Y, Z 

dimensionless coordinates 

m 


matrix (see eq. (16)) 

x,y,z 

a V a 2 

rectangular coordinates 
adjustable parameters 

j 

: \ 

V 

evaluation of quantity on 
fluid or wall side of 
contour 

@2 ’ ^3 

adjustable parameters 

Superscript: 


6 

variational operator 

T 


transpose of matrix 
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